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We perform a quantitative analysis of the nuclear modification factor in deuteron-gold collisions 
R dAu within the Color Glass Condensate approach, and compare our results with the recent data 
from RHIC experiments. Our model leads to Cronin enhancement at mid-rapidity, while at forward 
rapidities it predicts strong suppression of ft dAu at all pr due to low- a; evolution. We demonstrate 
that our results are consistent with the data for dAu charged hadron spectra, R dAu and R CP recently 
reported for rapidities in the interval r\ — ~ 3.2 by the BRAHMS experiment at RHIC. We also 
make a prediction for R pA at mid- rapidity in pA collisions at the LHC. 



PACS numbers: 24.85.+p, 12.38.-t, 12.38.Cy 

Recent observations QL H> El> S> 13 °f the suppression of high pt hadron yields at forward rapidities at RHIC have 
attracted considerable interest. The observed suppression is in sharp contradiction with the naive multiple scattering 
picture, in which the magnitude of Cronin enhancement observed at mid-rapidity is expected to increase further at 
forward rapidities, reflecting the growth of the number of scattering centers (partons) at small Bjorken x. On the 
t^J- . other hand, the observed effect has been predicted [f| 0, E as a signature of quantum evolution in the Color 



cttect has been predicted 

Glass Condensate (CGC) [l(J, EH E2j LCI LEI EH • Very recently, the first exploratory experimental results |ljj on 
the back-to-back azimuthal correlations of high pt particles separated by several units of rapidity in dAu collisions 
indicated the possible onset of the "mono-jet" behavior expected in the quantum CGC picture E3 (the azimuthal 
correlations in the classical approach to the CGC were studied in |l9|). 

Nevertheless, the origin of the observed effects is certainly not beyond a reasonable doubt at present; to clarify 
,-SZh ' it, one needs to perform dedicated and careful experimental and theoretical studies. While the data is qualitatively 
consistent with the predictions based on the CGC picture, a detailed comparison to the data requires a quantitative 
analysis taking into account, for example, the contributions of both valence quarks and gluons, and the influence 
D , of realistic fragmentation functions. Such an analysis is the goal of this note. Recently, related work in the more 
• traditional multiple scattering picture supplemented by shadowing has been done in Refs. [20ll2lT|. and in Ref. [2^ 
\ where the contribution of valence quarks scattering off the CGC has been addressed. 

In this work we use a simple model for the dipole-nucleus forward scattering amplitude which describes the onset of 
the gluon anomalous dimension in the Color Glass Condensate regime. Since the inclusive gluon and quark production 
cross sections in p{d)A collisions can be expressed in terms of the adjoint dipole-nucleus scattering amplitude, our 
model allows us to describe inclusive hadron production in deuteron-gold collisions at _J~s = 200 GeV at RHIC. Our 
model is based on a detailed analytical analysis performed in our previous publication T| stemming from the idea put 
forward in 

El- 
Inclusive cross section for production of a gluon in dA collisions was calculated in [2^ 0) E3 an d is given by 

da G A C F S A S d 1 



d?k dy a 8 tt (2tt) 3 fc 2 



d 2 z V 2 n G (z, Y-y) e~** V 2 N G (z, y), (1) 



where Sa and Sd are cross sectional areas of the gold and deuteron nuclei correspondingly and Y is the total rapidity 
interval. We assume a simple form of the scattering amplitude of the gluon dipole of transverse size zt = \z\ on the 
deuteron inspired by the two-gluon exchange 

na(z,lnl/x p ) = (1 - Xp) 4 x~ A 7ra 2 zfln(l/z T ^) — (2) 

Jd 

with A to be fixed later and x p the gluon's Bjorken x in the deuteron's (or proton's) wave function. Integrating over 
directions of z we rewrite Q as 



duff a s C F S A 



d 2 kdy 



f°° 1 

/ dz T J (k T z T ) In d ZT [z T d ZT N G (z T ,y)}, (3) 

Jo z tV 
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where \x is a scale associated with deuteron and is fixed at \x = 1 Gel/ thereof. The gluon djpole scattering amplitude 
on a gold nucleus Nq(zt, y) should be determined from the nonlinear evolution equation |lq|. Since an exact solution 
of the nonlinear evolution equation [lfij is a very difficult task we are going to construct a model for Na(zT,y) 
satisfying its asymptotic behavior: at Zt <C 1/Q s (y) one should have Ng(zt,v) ~ Zt, while at zt ^> 1/Q s (y) we 
should get N G (z T ,y) ~ 1 El HI 

is the nuclear saturation scale at rapidity y.) This behavior can be 

modeled by a simple Glauber-like formula 



N G (z T ,y) = 1 - exp 



-1(4^,4) 



(4) 



where 7(2/, zf.) will be given by (0. Note, that when 7 = 1 equations © and (@J reproduce the results of McLerran- 
Venugopalan model jl3L H4I |23| (for similar results see [2^ ) . 

At forward rapidities, in the deuteron fragmentation region, the Bjorken x of the nucleus acquires its lowest possible 
value for a given ^/~s, while the Bjorken x of the proton is close to unity. In that region rescatterings of valence quarks 
of the proton in a nucleus can give a substantial contribution to the hadron production cross section. This problem 
was discussed in a series of papers listed in (291 30] leading to the following expression for inclusive valence quark 
production cross section [3(| 



d< A S A 



d 2 k 2tt 



dz T z T J (k T z T ) [2- N Q (z T ,y)], (5) 



where Nq(zt,v) is the quark dipolc-nucleus forward scattering amplitude. In the quasi-classical approximation (7 
=1) Nq(zt, y) is given by the same quasi-classical formula with Q 2 (y) replaced by jf-Q 2 s (y) — ^Q 2 (y)- Therefore, 
by analogy with we model the quark dipole scattering amplitude Nq(zT,y) as 



NQ(z T ,y) = 1 - exp 



AV T N r Qs 



(6) 



To model the anomalous dimension j(y, zip) we use the following interpolating formula 




ay, 



2 • 



t(»,4) = ; 1 1+ efa , 4)+ ;^ 4)+7C(3)c i . < 7 > 



where 



(V 2 )(y-yo) 

and c is a constant to be fitted. This form of the anomalous dimension is inspired by the analytical solutions to the 
BFKL equation [3]]]. Namely, in the limit zt — > with y fixed we recover the anomalous dimension in the double 
logarithmic approximation 7 w 1 — yT/(2£). In another limit of large y with zt fixed, Eq. (J7J reduces to the expression 
of the anomalous dimension near the saddle point in the leading logarithmic approximation 7kU ^, . . Therefore 
Eq. Q mimics the onset of the geometric scaling region |23,|23J. A characteristic value of zt is zt ~ 1/(2 hr), so we 
will put 7(y,4) «7(y. 1/(4 *!)). 

The saturation scale Q s (y) that we use is the same as the one used in [34| to fit the low-x DIS data and in [35j to 
describe the hadron multiplicities at RHIC. It is given by 

Q 2 s (y) = A 2 A 1 / 3 e Xy = 0.13 GeV 2 e A y N coU . (9) 

Here N co n is the number of binary collisions at a given centrality in a dAu collision. Parameters A = 0.6 GeV and 
A = 0.3 are fixed by DIS data |3J|. The initial saturation scale used in JHJ is defined by Q 2 s0 — Q 2 {yo) with yo the 
lowest value of rapidity at which the low-x quantum evolution effects are essential. 

The Cronin effect [2^ is usually attributed to multiple rescatterings of partons in the nucleus 00, However, 
it is also present in the low energy data, i.e. at energies where saturation is unlikely to play a significant role for the 
production of high pt particles. For example, at y/s = 20 GeV the nuclear enhancement for ir^ produced in proton- 
nucleus collisions peaks at kx — 4 GeV |3fij | . This implies that the typical nonperturbative scale k associated with such 
low energy hadronic rescatterings may be rather large. It becomes much smaller than Q s (y) at high energies/rapidities 
as one can see from However, at the central rapidity region at RHIC the influence of this non-perturbative scale 
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cannot yet be neglected. To take it into account in describing the nuclear modification at RHIC we shift the saturation 
scale in the Glauber exponents and © as follows Q 2 — ► Q 2 + n 2 A 1 / 3 . This shift is also performed for Q 2 in 
(jHJ- In our numerical calculation we chose two values of n: k = 1 GeV takes into account additional momentum 
broadening due to a nonperturbative effects and n — neglects such effects. 
The nuclear modification factor is usually defined as 

dN dAu 



RdAu(k T ,y)= AT d2 k dNp ~ > ( 10 ) 

iVco11 d 2 kdy 

where d j% k dy and Ji^^ are multiplicities of hadrons per unit of phase space in dAu and pp collisions. Both expressions 
for gluon © and quark J^} production contribute to the hadron production cross section in dAu collisions. The cross 
section of hadron production is given by 

dadGA -(k T /z)Df rag (z,k T )F(k T /z,y) 



cPkdy J z 2 d 2 kdy 

f dz daf) A n 
+ J ~z* ~&k {kT,z) Xqv{y ' kT,z) D i^ { '^ kr) F{kr/Z > V) - 

We use the LO fragmentation functions from Ref. [23] • We choose the renormalization scale of the fragmentation 
functions to be kx- Equation JSJ is derived for production of a valence quark in the deuteron fragmentation region. 
To generalize it to smaller values of Bjorken x one has to convolute it with the deuteron's valence quark distribution, 
which is fixed by quark counting rules at high x and by the leading Regge trajectory at low x 

xq v (x) = 1.09(l-x p ) 3 x° p 5 , (12) 

where x p = (h? / \/~s) e n . (Eq. (|12(l is normalized to give the distribution of a single valence quark in the deuteron 
to keep normalization the same as in Q.) Valence quarks are increasingly less important at low x [39l |. where the 
quark production is dominated by gluons splitting in qq pairs. The factor of x^ 5 insures that this is indeed the case 
here |3!j • Analogously, the high x behavior of the nuclear gluon distribution is taken into account by introducing the 
function F{kx, y) 

( A 2 \ 13a ° 

F(k T ,y) = (l-x A r fp^J • (13) 

where the Bjorken x of a gluon in the nuclear wave function is given by xa — {kr / ' \fs) and a s = 0.3. The 
last factor in Eq. (|13[1 arises when we impose momentum conservation constraint on the anomalous dimension of the 
distribution functions. Namely, we use the following phcnomcnological parametrization of the anomalous dimension 
in the Mellin momentum variable u> |32j| 

7 (w) =a 8 (^-lj. (14) 

This parametrization takes into account high x corrections to the QCD splitting functions. 

The differential hadron multiplicity can be calculated by dividing Ijllll by the total inelastic cross section <7dAu l° r 
a given centrality selection. The baseline pp multiplicity is calculated by expanding the Glauber exponent (QJ to the 
leading term at Zt <C 1/Q s - The free parameters of our model are yo in ©, which sets the initial value of y at 
which the quantum evolution sets in, c in J7J), which describes the onset of quantum regime, the momentum scale k, 
which specifies the typical hadronic rescatterings momentum, and /i in which is the infrared cutoff. The value of 
fj, = 1 GeV and the range of values for k = -r 1 GeV are fixed by lower energy data. Parameters yo and c are fitted 
to RHIC dAu data reported by BRAHMS collaboration [ij, |2( . The parameter A from 10 is fixed by the DIS data 
and is not a free parameter in our model. 

The data reported in Ref. 0] is for charged particles at pseudo-rapidities r/ = 0, 1 and for negative ones at pseudo- 
rapidities r) — 2.2,3.2. At forward rapidities (in the deuteron fragmentation region) the valence quarks begin to 
dominate over gluons in the production of hadrons with high transverse momenta. In particular, in pp collisions 
this leads to an asymmetry between positive and negative hadrons - an effect which is well-established (see ^(J an( l 
references therein). Since the nuclear modification factor RdA has been experimentally defined as the ratio of dAu 
and pp cross sections, this factor is modified by the isospin asymmetry effects. Unfortunately, it is difficult to evaluate 
quantitatively the magnitude of these effects - the isospin dependence of fragmentation functions is poorly known, 
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and the relative importance of valence quarks and gluons in various kinematical regions heavily depends on the choice 
of the structure functions. Nevertheless, to account for the influence of this effect we performed calculations for two 
limiting cases: (i) assuming no isospin dependence for the valence quark fragmentation and (ii) in the opposite limit 
of the constituent quark model, with u-quarks fragmenting only into positive hadrons and d-quarks fragmenting only 
into negative ones. 

The results of our calculations are presented in Figs. ^ HI and 01 along with the data collected by BRAHMS 
collaboration 0,0. In these figures we use c = 4 with y — 0.6 for both k = and k = 1 GeV. We would like to 
emphasize that the ratios RdAu and Rcp are almost insensitive to the values of k and fj, at 77 > 1 and pt > 2 GeV. 
Their dependence on yo is also weak at forward rapidities. 

In Fig. n we present our calculation of the charged particle transverse momentum spectra in dAu collisions at 
several different rapidities compared to BRAHMS data 2] . We find a reasonable agreement with experimental data 
Q. To evaluate the degree of agreement with the data one should also keep in mind that BRAHMS data at r\ = 
and 1 are for the average of positive and negative hadrons, and so contain the baryons which production dynamics 
still remains puzzling at present (the shown data at 77 = 2.2 and 3.2 are for the negative hadrons only). In Fig. El 
we show the nuclear modification factor R dAu as a function of pr at different rapidities calculated in our model and 
compared to the data from |2J . At rapidity r/ < ?/o w e observe Cronin enhancement of the nuclear modification factor 
at pt ~ 2 3 GeV due to the multiple rescatterings of the deuteron in the gold nucleus El IH, |37} • At 77 > yo 
the low- a; quantum evolution effects modify the anomalous dimension 7 leading to suppression in R dAu at forward 
rapidities and disappearance of the Cronin maximum in accordance with our qualitative predictions in Ref. 7] (see 
also 0]). Fig. demonstrates the centrality dependence of the hadronic spectra by plotting the central-to-peripheral 
ratio R for the same rapidities as in Figs. ^ and [21 It is important to emphasize that Rcp is much less sensitive 
than RdAu to isospin-dependent effects. At mid-rapidity the Cronin maximum increases and moves to the right as 
centrality increases. Conversely, at forward rapidities, r\ 3> yo, the suppression gets stronger with centrality since, 

approximately, R dAu ~ yJl/N^t 00. This behavior at forward rapidities is in agreement with the BRAHMS data 

S B. Nuclear modification factor obtained by numerical solution of the BK equation 0] as well as other approaches 
|30| qualitatively agree with our conclusions. Further research in the area included an analysis of running coupling 
corrections fill ] and a study of similar suppression in di-lepton production |42j| . 

It will be interesting to check the predictions of our approach at the LHC energy of ^fs = 5.5 TeV. In Fig. 01 we 
show our result for the nuclear modification factor R pA for pA collisions at LHC at mid-rapidity, compared to R dA 
for RHIC dAu collisions at 77 = 3.2. Our model predicts that the nuclear modification factor will be quite similar 
for both cases. This conclusion seems natural to us since the effective values of nuclear Bjorken x for mid-rapidity 
LHC collisions will be similar to the effective x achieved in the forward rapidity RHIC collisions. If observed, the 
mid-rapidity suppression predicted in Fig. 01 for pA collisions at LHC would indicate that an even stronger suppression 
due to the CGC initial state dynamics should be present in AA collisions at LHC. The overall high-p^ suppression in 
mid-rapidity R AA at the LHC would then be due to both the initial state saturation/CGC dynamics and jet quenching 
in quark-gluon plasma |43l E EH l46l | . 

In summary, we presented a simple but quantitative model which incorporates the main features of the small- a; 
evolution in the Color Glass Condensate for particle production in ultra relativistic proton (deuteron)-heavy ion 
collisions. We find that at ^/s = 200 GeV the evolution sets in at rapidity yo ~ 0.2. As a result, at central rapidities 
the nuclear modification factor R dAu exhibits Cronin enhancement at kx — 2 3 GeV, whereas at forward rapidities 
it is strongly suppressed at all pt ■ 
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